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ABSTRACT 

We investigate how the fraction of broad-line sources in the AGN population changes with X-ray 
luminosity and redshift. We first construct the rest-frame hard-energy (2 — 8 keV) X-ray luminosity 
function (HXLF) at z = 0.1 — 1 using Chandra Lockman Hole-Northwest wide-area data, Chandra 
Deep Field-North 2 Ms data, other Chandra deep field data, and the ASCA Large Sky Survey data. 
We find that broad-line AGNs dominate above ~ 3 x 10 43 ergs s _1 and have a mean luminosity of 
1.3 x 10 44 ergs s _1 . Type II AGNs can only become an important component of the X-ray population 
at Seyfert-like X-ray luminosities. We then construct z = 0.1 — 0.5 and z — 0.5 — 1 HXLFs and compare 
them with both the local HXLF measured from HEAO-1 A2 survey data and the z = 1.5 — 3 HXLF 
measured from soft-energy (0.5 — 2 keV) Chandra and ROSAT data. We find that the number density of 
L x > 10 44 ergs s _1 sources (quasars) steadily declines with decreasing redshift, while the number density 
of L x = 10 43 — 10 44 ergs s _1 sources peaks at z = 0.5 — 1. Strikingly, however, the number density of 
broad-line AGNs remains roughly constant with redshift while their average luminosities decline at the 
lower redshifts, showing another example of cosmic downsizing. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: formation — galaxies: 
active 



1. INTRODUCTION 

It is becoming increasingly clear that low-redshift accre- 
tion may be the dominant phase of supermassive black hole 
(SMBH) growth (Barger et al. 2001; Cowie et al. 2003). 
Although theoretical models of SMBH formation have re- 
lied on comparisons with the optical quasar luminosity 
function, optically accessible broad-line active galactic nu- 
clei (AGNs) form only about 30% (Barger et al. 2003b) of 
the X-ray background (XRB). Explanations for the lack of 
broad lines in most AGN optical spectra include absorp- 
tion along the line-of-sight or dilution of the emission-line 
signatures by stellar light (Moran, Filippenko, & Chornock 

2002) . Regardless of the reason, optically and soft X-ray 
selected samples substantially undercount the AGN pop- 
ulation relative to hard X-ray selected samples, at least at 
intermediate luminosities (Cowie et al. 2003; Barger et al. 
2003a). 

The hardest band for which the hard X-ray luminosity 
function (HXLF) can presently be determined is rest-frame 
2 — 8 keV. The » 2 Ms exposure of the Chandra Deep 
Field-North (CDF-N) samples a large, distant cosmolog- 
ical volume down to a very faint hard X-ray flux limit 
(/2-8 koV ~ 1-4 x 10~ 16 ergs cm~ 2 s -1 ; Alexander et al. 

2003) . However, this is substantially fainter than the X- 
ray flux of the sources that contribute the most to the XRB 
(Cowie et al. 2002), and the low-redshift volume sampled 
is small. To sample a large cosmological volume at low red- 
shifts, we have mapped with Chandra an ~ 0.4 deg 2 area in 
the Lockman Hole- Northwest (LH-NW) field to a flux level 
of f 2 -s kcv ~ 3 x 10~ 15 ergs cur 2 s" 1 (Yang et al. 2003), 
where most of the hard XRB is resolved. In this paper we 



use our LH-NW data, together with deep Chandra data 
on the CDF-N, A370, SSA13, and SSA22 fields and bright 
ASCA data from Akiyama et al. (2000), to construct rest- 
frame 2 — 8 keV luminosity functions at low redshifts. We 
then investigate the evolution with redshift of AGN pop- 
ulations separated by optical spectral type. We assume 
Q M = 1/3, fi A = 2/3, and H Q = 65 km s" 1 Mpc" 1 . We 
use L x to denote the rest-frame 2 — 8 keV luminosity. 

2. X-RAY SAMPLE SELECTION 

The LH-NW data are taken from the X-ray and optical 
catalogs of Y. Yang et al. (in preparation) and A. T. Stef- 
fen et al. (in preparation). Here we consider only sources 
with fluxes above 5 x 10 -16 ergs cm -2 s _1 (0.5 — 2 keV) and 
5 x 10~ 15 (2 — 8 keV) to provide good, high-significance 
samples. We also only consider those sources that have 
been targeted spectroscopically. The only source property 
we considered when making our masks for spectroscopy 
was X-ray flux (sources with higher X-ray fluxes had higher 
priority when we had to choose one X-ray source over an- 
other due to mask conflicts). We spectroscopically ob- 
served 150 sources out of our total 2 — 8 keV sample of 
228. We were able to confidently identify 100 of these 
150 sources (67%) based on multiple emission and/or ab- 
sorption line features. Two of the identified sources are 
stars. A practical optical limit for obtaining redshifts 
for X-ray sources using the Deep Extragalactic Imaging 
Multi-Object Spectrograph (DEIMOS; Faber et al. 2002) 
on Keck II is about R = 24 (see Barger et al. 2003b). 
Other spectra may not have been identifiable due to slits 
falling too close to an edge of the CCD or bright neighbor 
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Fig. 1. — Redshift (solid for spectroscopic, open for photometric) vs. 2 — 8 kcV flux for the X-ray samples (squares — CDF-N; triangles — 
A370, SSA13, and SSA22; diamonds — LH-NW; circles — ASCA). Broad-line optical spectra are denoted by large solid symbols. Unidentified 
sources are denoted by small solid symbols at z < 0. Solid curves show flux vs. redshift for L x = 10 ergs s" - 1 (lowest curve), 10 43 , 10 44 , 
and 10 45 (highest curve), computed with a K-correction for a T = 1.8 power law spectrum (see Barger et al. 2002). 



objects compromising the spectra. 

The faintest X-ray data are taken from the w 2 Ms CDF- 
N catalogs of Alexander et al. (2003; X-ray) and Barger et 
al. (2003b; optical). Excluding stars and the small number 
of extended X-ray sources, the 2 — 8 keV sample contains 
326 sources, 68% of which have either spectroscopic (165) 
or reliable photometric (56) rcdshifts. Our other Chandra 
and ASCA samples are summarized in Cowie et al. (2003). 
Hereafter, we refer to the A370, SSA13, and SSA22 Chan- 
dra samples collectively as "other" . 

In our analysis we place every X-ray source that was not 
spectroscopically identified for any reason into each and 
every redshift interval in order to obtain a strict upper 
bound on the number of AGNs that could lie in that red- 
shift interval, modulo the existence of a substantial pop- 
ulation of Compton-thick sources not detected by Chan- 
dra. We note in passing that Barger et al. (2002) found 
from the 1 Ms CDF-N data that a substantial popula- 
tion of Compton-thick sources is not needed to explain 
the greater than 8 keV XRB, provided that the ASCA 
and BeppoSAX XRB measurements, instead of the lower 
HE A O-l A2 measurement, are assumed at energies less 
than 8 kcV. For the LH-NW field, we only ever consider 
the sources that were spectroscopically observed (identi- 
fied or not). This is acceptable because the spectroscopy 
are purely from targeted observations of the X-ray sam- 
ples. For the other fields, the spectroscopy are not purely 
from targeted observations of the X-ray samples, so we 
need to include all of the sources. 

We determined the solid angle covered by the com- 
bined samples at a given flux by comparing the observed 



numbers of sources versus flux with the averaged number 
counts in the 2 — 8 keV band from Cowie et al. (2002). 
This method allows a simple treatment of the incomplete- 
ness that was modeled in computing the counts. However, 
the counts in Cowie et al. also include the low CDF-S 
counts, which may affect the normalization at the 10% 
level relative to the LH-NW average (Yang et al. 2003). 
We consider this to be a reasonable estimate of the sys- 
tematic errors in the present analysis. The solid angle 
covered by the combined 2 — 8 keV samples ranges from 
0.002 deg 2 at the faintest fluxes to 5.8 deg 2 at the highest 
fluxes. At 2.3 x 10~ 14 ergs cm~ 2 s" 1 (2-8 keV) the solid 
angle is 0.52 deg 2 . 

In Figure 1 we show redshift versus 2 — 8 keV flux for the 
ASCA (circles), LH-NW (diamonds), CDF-N (squares), 
and other (triangles) samples. The LH-NW sources nicely 
fill in the almost an order of magnitude flux gap between 
the CDF-N and ASCA samples. The solid curves corre- 
spond to loci of constant L x . Any source more luminous 
than L x = 10 42 ergs s _1 is very likely to be an AGN on 
energetic grounds (Zezas, Georgantopoulos, & Ward 1998; 
Moran et al. 1999), though many of the intermediate lu- 
minosity sources do not show obvious AGN signatures in 
their optical spectra. 

3. THE AGN X-RAY LUMINOSITY FUNCTION 

Cowie et al. (2003) used Chandra, ASCA, and ROSAT 
data to determine rest-frame 2 — 8 keV AGN luminosity 
functions in two redshift intervals, z = 0.1 — 1 and z = 2—4. 
They did not have very many sources in their sample with 
2-8 kcV fluxes - 2 x 10~ 14 to 10~ 13 ergs cur 2 s" 1 since 
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that intermediate flux range is not well covered by either 
small-area, deep Chandra observations or bright ASCA ob- 
servations. With our wide-area LH-NW data, we now have 
good coverage of the intermediate flux range and can there- 
fore improve the determination of the z = 0.1 — 1 HXLF. 
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Fig. 2. — (a) Rest-frame 2 — 8 keV luminosity function per unit 
logarithmic luminosity for z = 0.1 — 1 (solid diamonds), computed 
from observed-frame 2 — 8 keV assuming an intrinsic T = 1.8, for 
which there are only small differential ii"-corrections to rest-frame 
2 — 8 keV. No correction was made for instrinsic absorption, but 
in rest-frame 2 — 8 keV, average absorption corrections are small. 
Maximal z = 0.1 — 1 HXLF (open diamonds) was found by assign- 
ing all of the unidentified sources redshifts at the volume-weighted 
centers of each and every rcdshift interval. True HXLF lies some- 
where between the solid and open diamonds. Open squares show 
the z = 0.1 — 1 HXLF for broad-line AGNs only. Poissonian la 
uncertainties are based on the number of sources in each logarith- 
mic luminosity bin. (b) Rest-frame HXLF from (a) broken into two 
smaller rcdshift intervals, z = 0.1—0.5 (solid squares) and z = 0.5—1 
(solid diamonds) . Maximal HXLFs are shown with open symbols. 
The 2 = 1.5 — 3 HXLF (dashed line; computed from observed-frame 
0.5 — 2 keV) and the local Seyfert galaxy HXLF (solid line; from 
Piccinotti et al. 1982) are shown for comparison. Dotted line is 
maximal z = 1.5 — 3 HXLF, found by assigning redshifts to all of 
the unidentified sources at the volume- weighted centers of each and 
every redshift interval. True z = 1.5 — 3 HXLF lies somewhere 
between the dashed and dotted lines. 

In Figure 2a we show our measured HXLF for the 
z = 0.1 — 1 interval (solid diamonds), computed following 
Cowie et al. (2003), who used the traditional 1/V a method 
of Felten (1977). The uncertainties are Poissonian, based 
on the number of galaxies in each logarithmic luminosity 



bin. However, since incompleteness is a potentially larger 
source of error, we also recomputed the HXLF (open dia- 
monds) by assigning all of the unidentified sources to the 
volume-weighted centers of each and every redshift inter- 
val. Because the spectroscopic identifications arc much 
more complete at the high X-ray fluxes, the associated 
systematic uncertainties are larger at low L x . The HXLF 
of only the broad-line AGNs is shown by the large open 
squares in Figure 2a. In the high luminosity range where 
the HXLF is steep, we can see that the broad-line AGNs 
dominate the total HXLF. 

The uncorrected z = 0.1 — 1 HXLF can be well repre- 
sented by a conventional broken power law fit of the form 
6.6 x 10~ 6 L44 137 per Mpc 3 above 2 x 10 43 ergs s" 1 and 
1.5 x 10" 5 L^ - 87 below, where £44 is the luminosity in 
units of 10 44 ergs s _1 . The integrated z = 0.1 — 1 light 
density is 1.8 x 10 39 ergs s" 1 Mpc~ 3 . Half of the total un- 
corrected light in the z = 0.1 — 1 interval arises in sources 
more luminous than 2.1 x 10 ergs s , and about one- 
third arises in broad-line AGNs. 

In Figure 2b we show our measured HXLF over the lu- 
minosity range L x = 10 43 to 10 45 ergs s _1 , where most of 
the sources have been spectroscopically identified, divided 
into two smaller redshift intervals, z = 0.1 — 0.5 (squares) 
and z = 0.5 — 1 (diamonds). For comparison, we show the 
local Seyfert galaxy HXLF (solid line) from Piccinotti et 
al. (1982; converted to H = 65 km s _1 Mpc -1 and, using 
their assumed T = 1.7, to 2 — 8 keV) and the high-redshift 
HXLF (dashed line) from Cowie et al. (2003; updated to 
include our LH-NW data and the 2 Ms CDF-N data, and 
to cover the interval z = 1.5 — 3). We have computed the 
maximal high-redshift HXLF (dotted line) by assigning all 
of the unidentified sources to the volume-weighted cen- 
ters of each and every redshift interval. The high-redshift 
HXLF is much tighter than in Cowie et al. (2003) because 
of our slightly lower rcdshift interval and our higher spec- 
troscopic completeness. From Figure 2b we see that the 
number density of L x > 10 ergs s _1 sources steadily 
declines with decreasing rcdshift. However, the number 
density of L x = 10 43 — 10 44 ergs s _1 sources is largest in 
the z — 0.5 — 1 interval and smaller at both higher and 
lower redshifts (see also Fig. 4). 



4. X-RAY FLUXES AND LUMINOSITIES 

Prior to Chandra and XMM-Newton, various groups 
(e.g., Comastri et al. 1995; Gilli, Salvati, & Hasinger 2001) 
using XRB population synthesis models based on unified 
AGN schemes predicted that a mixture of unobscured type 
I AGNs and intrinsically obscured type II AGNs were 
needed to produce the bulk of the XRB. The obscured sys- 
tems were expected to lie at high redshifts (z = 2 — 3), with 
the high luminosity end (type II quasars) contributing sub- 
stantially to the XRB. Contrary to the predictions, optical 
observations of the sources detected in deep Chandra and 
XMM-Newton hard X-ray surveys have found a redshift 
distribution that peaks at z < 1 (Barger et al. 2002, 2003b; 
Hasinger 2003). Only a handful of type II quasars (here 
we define quasars as having L x > 10 44 ergs s _1 ) have been 
detected, and these contribute only a small fraction of the 
XRB (e.g., Barger et al. 2003b). 

In Figure 3 we plot for the redshift intervals z = 0.1 — 1 
and z = 1.5 — 3 the fractional number of X-ray sources in 
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each luminosity bin, by optical spectral type, versus rest- 
frame 2 — 8 kcV luminosity. From Figure 3a we see that 
the dominant population at quasar luminosities, where 
our redshift identifications are very complete, is broad- 
line AGNs. However, at more Seyfcrt-likc luminosities 
(L x < 10 44 ergs s _1 ), the broad- line AGN fraction de- 
clines rapidly with decreasing luminosity. Figure 3b shows 
a similar trend with luminosity. These results suggest a lu- 
minosity dependence in optical spectral type: type I AGNs 
dominate at quasar luminosities, while type II AGNs can 
only become an important component of the X-ray popu- 
lation at Seyfert-like luminosities. 
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Fig. 3. — Fractional number of X-ray sources by optical spec- 
tral type (dark shading — non-broad-line sources; medium shading — 
broad-line AGNs; light shading — unidentified sources) versus 2 — 
8 kcV luminosity for the redshift intervals (a) z = 0.1 — 1 and (b) 
z = 1.5 — 3. Luminosities were computed with a if-correction for 
a r = 1.8 power law spectrum. All of the unidentified sources were 
assigned rcdshifts at the volume-weighted centers of each and every 
redshift interval. 



To quantify this, we note that in the z = 0.1 — 1 redshift 
interval there are 20 sources with quasar luminosities, of 
which 15 are broad-line AGNs. Only three of the uniden- 
tified sources would have quasar luminosities if placed in 
this redshift interval, so the systematic errors are small, 
and the bulk of the quasars are broad-line AGNs. In the 
z = 1.5 — 3 interval there are 65 sources with quasar lumi- 
nosities, of which 58 are broad-line AGNs. There are 56 
unidentified sources that would have quasar luminosities 



if placed in this redshift interval, so the systematic uncer- 
tainties are larger. However, at least half (and probably 
considerably more) of the sources with quasar luminosities 
in this redshift interval are broad-line AGNs. 

5. DISCUSSION 

Barger et al. (2003a) graphically showed how the L x = 
10 43 — 10 44 ergs s _1 (Seyfert-like luminosities) number den- 
sity slowly increases with decreasing redshift from z — 6 to 

-i i) i » s 1 (quasar imiii 
1.5 - 3. They also 



z = 0.1, while the L x = 10 44 - 10"" ergs s 



nositics) number density peaks at z 
showed that the evolution of the L x = 10 44 — 10 45 ergs s _1 
number density closely matches in shape the evolution of 
optically selected quasars. Given these trends and our 
interesting result from Figure 3 that while the dominant 
population at quasar luminosities is broad-line AGNs, at 
lower luminosities the broad-line fraction drops dramat- 
ically, here we explore a different aspect of the number 
density evolution. 

In Figure 4 we show number density versus redshift for 
the sources with L x > 10 43 ergs s _1 . Rather than di- 
viding the sources into luminosity classes, as was done in 
Barger et al. (2003a), we have instead divided the sources 
into broad-line (open squares) and non-broad-line (solid 
circles) AGNs. The horizontal bars show the maximal 
number densities for the non-broad- line AGNs, obtained 
by assigning all of the unidentified sources redshifts at the 
volume- weighted centers of each and every redshift bin and 
then retaining only those with L x > 10 43 ergs s _1 . The 
bars are not consistent with one another because all of the 
unidentified sources are included in each and every redshift 
bin, as long as their luminosities at those redshifts exceed 
10 43 ergs s _1 
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Fig. 4. — Number density of sources with L x > 10 4a ergs s _1 
vs. redshift. Broad-linc (non-broad-line) AGNs are denoted by 
open (solid) symbols. Points above (below) z = 2 were determined 
from observed-frame 2 — 8 keV (0.5 — 2 kcV), assuming an intrin- 
sic r = 1.8. Poissonian la uncertainties are based on the number 
of sources in each redshift bin. Horizontal bars show the maximal 
number densities for the non-broad-line AGNs, found by assigning 
redshifts to all of the unidentified sources at the volume-weighted 
centers of each and every redshift bin and then retaining only those 
with L x > 10 43 ergs s — 1 . Dashed line shows the mean luminosity 
(right-hand axis) of the broad-line AGNs with redshift. 

From Figure 4 we see that the number density of non- 
broad-line AGNs with L x > 10 43 ergs s _1 rises with dc- 
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creasing redshift, until the last rcdshift bin. Strikingly, 
however, the number density of broad-line AGNs remains 
roughly constant with rcdshift while their average lumi- 
nosities (dashed line; right-hand axis) decline at the lower 
rcdshifts. This is another example of cosmic downsizing. 
In the terminology of previous quasar studies, the broad- 
line AGNs are undergoing pure luminosity evolution, not 
density evolution. Pure luminosity evolution is consis- 
tent with the most recent optical determinations (Boyle 
et al. 2000). The increase in the number density of non- 
broad-line sources and the decrease in the broad-line AGN 
average luminosities results in non-broad-line sources in- 
creasingly dominating the light at lower redshifts, modulo 



the existence of a substantial population of Compton-thick 
sources. 
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